Long-range gravitational-like interaction in a neutral atomic cold gas 
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A quasi resonant laser induces a long-range attractive force within a cloud of cold atoms. We take 
advantage of this force to build in the lab quasi ID systems of particles with a gravitational-like 
interaction, at a fluid level of modeling. We give experimental evidences of such an interaction in a 
Strontium cold gas, studying the density profile of the cloud, its size as a function of the number of 
atoms and its breathing oscillations. 
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When interactions between the microscopic compo- 
nents of a system act on a length scale comparable to 
the size of the system, one may call them "long range" . 
For instance, the inverse-square law of the gravitational 
force between two point masses is one of the most cel- 
ebrated and oldest laws in Physics. In the many parti- 
cles world, it is responsible for dramatic collective effects 
such as the collapse of a gravothermal catastrophe [ij 
or the gravitational clustering which is the main mecha- 
nism leading to the formation of the structure of galax- 
ies in the present universe. Beyond gravitation, such 
long range interactions are present in various physical 
fields, either as fundamental or as effective interactions: 
plasma physics, two dimensional fluid dynamics, wave- 
particles systems. . . They deeply influence the dynami- 
cal and thermodynamical properties of such systems. At 
the thermodynamic equilibrium, long range interaction 
is at the origin of very peculiar properties, especially for 
attractive systems: the specific heat may be negative; 
canonical (fixed temperature) and microcanonical (fixed 
energy) ensembles are not equivalent. These special fea- 
tures have been known for a long time in the astrophysics 
community, in the context of self gravitating systems. 

After the seminal works of Lynden-Bell and Wood Q 
and Thirring Q, many contributions followed on this 
subject (see for instance (3] for a recent review), so that 
the equilibrium characteristics of attractive long range in- 
teracting systems are theoretically well established. This 
situation is in striking contrast with the experimental 
side of the problem: there is currently no controllable 
experimental system exhibiting the predicted peculiar- 
ities. There have been some proposals to remedy to 
this situation: O'Dell et aL [5[ have suggested to cre- 
ate an effective 1/r potential between atoms in a Bose- 
Einstein condensate using off-resonant laser beams; re- 
cently, Dominguez et al. have proposed to take advantage 
of the capillary interactions between colloids to mimic 
two-dimensional gravity 0, and Golestanian has sug 
gested experiments using thermally driven colloids 



However, these proposals have not been implemented yet, 
and so far the dream of a tabletop galaxy remains elusive. 
The key results of this paper are to show some experi- 
mental evidences of a gravitational-like interaction in a 
quasi one-dimensional (hereafter ID) test system consist- 
ing in a cold gas of Strontium atoms in interaction with 
two counter-propagating quasi-resonant lasers. 

First let us show that the quasi ID {cold atomic gas 
-|- ID quasi-resonant laser beams} system can indeed 
exhibit an ID gravitational-like force. For that pur- 
pose, we consider an atomic gas with a linear density 
n{z) in interaction with two counter-propagating laser 
beams. The two beams intensities I+{z) and I-{z), 
where /+(— oo) = /_(-|-oo) = Jo, respectively propa- 
gating in the positive and negative direction, are much 
smaller than the atomic line saturation intensity Is ■ Thus 
the atomic dipolar response is linear. The radiation pres- 
sure force of the lasers on a single atom having a longi- 
tudinal velocity Vz is given by [8|: 

F±{z,vz)^±hk^——^^^^^—^^^, (1) 



2 A{STkvz)^ + T'^ Is 

where h is the reduced Planck constant, F the bare 
linewidth of the atomic transition, k the wavenumber 
and 5 the frequency detuning between an atom at rest 
and the lasers. The attenuation of the laser intensity is 
given by: 

d^± = T:^n(^)dz, (2) 
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is the absorption cross-section for a cloud of atoms having 
a normalized longitudinal velocity distribution f(i'z), and 
27rL^ is the transverse section of the cloud. The optical 
depth is defined as 
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n{z)dz . 
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Atoms also experience a velocity diffusion due to the 
random photon absorptions and spontaneous emissions. 
In experiments, S < such that the force, given in 
Eq. dD), is a cooling force counteracting the velocity dif- 
fusion. In this context, we can derive a fluid equation 
starting from the appropriate ID Fokker-Planck equation 
for the one-point probability density function f(z,Vz,t) 
describing particles in a trap with a longitudinal trap 
frequency uJz and interacting with the pair force ^ . We 
assume that / may be written as 
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u{z, t) 



(5) 



where f{vz) is even and homogeneous along the sample; 
M is a (small) macroscopic velocity and the velocity dis- 
persion is characterized by a time modulation A{t). In- 
tegrating the Fokker-Planck equation over Vz dvz and 
looking for a stationary solution we obtain: 
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sgn(s — z)n{s)ds = . (6) 



m the atom mass, ks the Boltzmann constant and Vz 
the rms velocity. The temperature, defined as fc^T = 
mVz = m Vzf{vz)dvz, results from the equilibrium 
of cooling forces and the velocity diffusion. It appears in 
the first term of Eq. (|6]) which governs the spreading of 
the density distribution. One can note that, in contrast 
with the 2D and 3D cases, this term always prevents 
the collapse of the cloud. The second term is due to 
the external holding trap, supposed here to be harmonic. 
The long range attractive interaction appears in the third 
term. Eq. (j6]) is derived in the weak absorption limit, i.e. 
assuming 6^1. One defines a characteristic interaction 
length as 



keT 1 
Cm N 



(7) 



with N is the number of atoms, C = -^^^r^if, and rj = 



+ 00 



oo r^+4(a-fct»,)^ K^^)'^"^ • The self-gravitating regime is 
then reached if Li <C Lni where 
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(8) 



trap 



is the characteristic length of the non-interacting gas in 
its external harmonic holding potential. The quantities 
ZR and kBTtrap/2 are respectively the Rayleigh length 
and the depth of the trapping potential. If the inequality 
Li <^ Lni is fulfilled, Eq. ([6]) yields the profile of a ID 
self-gravitating gas in thermal equilibrium 
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The attractive force coming from the beams absorp- 
tion, as described above, is known since the early days of 



laser cooling and trapping [lO[. Unfortunately, in a 3D 
setting this attractive force is usually dominated by the 
repulsive force due to photons reabsorption [llj, which 
in the small b limit may be seen as an effective repulsive 
Coulomb force. By contrast, in a ID configuration with 
an elongated cloud along the cooling laser beams, the 
probability of photons reabsorption is reduced by a factor 
of the order of L^/L^, in comparison with the isotropic 
cloud having the same longitudinal optical depth. 
and Lz are respectively the rms transverse and longi- 
tudinal size of the cold cloud. In our experiment, the 
reduction factor is about 2 x 10"^, so that the repulsive 
force can be safely ignored. 

The sample preparation is done in the same way as de- 
picted in [12| : more details about laser cooling of Stron- 
tium in a magneto-optical trap (MOT) can be found in 
[isj . After laser cooling, around 10^ atoms at T ~ 3 /iK 
are loaded into a far detuned dipole trap made of a 
120 mW single focused laser beam at 780 nm. We directly 
measure the longitudinal trap frequency uJz = 6.7(0.5) Hz 
from relaxation oscillations of the cold cloud (see exam- 
ple of temporal evolutions in Fig. [T]). The radial trap 
frequency a;j_ = 470(80)IIz is deduced from cloud size 
measurements. The beam waist is estimated at 23(2) jjxa. 
leading to a potential depth of Ttiap — 20 /^K. 

50 ms after loading the dipole trap (corresponding to 
t = in Fig. [TJ, a counter- propagating pair of laser 
beams, red-detuned with respect to the ^5*0 ^Pi in- 
tercombination line at 689 nm (radiative lifetime: 21 /xs), 
are turned on for 400 ms. These beams, aligned with re- 
spect to the longitudinal axis of the cloud, generate the 
effective ID attractive interaction. When the ID lasers 
are on, we apply a B — 0.3 G magnetic bias field, for two 
important reasons: First, the Zeeman degeneracy of the 
excited state is lifted such that the lasers interact only 
with a two-level system made out of the m = — >■ m = 
transition which is unsensitive to the residual magnetic 
field fiuctuation. Second, the orientation of magnetic 
field bias, with respect to the linear polarization of the 
dipole trap beam, is tuned to cancel the clock (or tran- 
sition) shift induced by the dipole trap on the transition 
of interest [l^ . 

The temperature along the ID laser beams, in our ex- 
perimental runs, is found to be in the range of 1 — 3^K. 
Even at such low temperatures, and in sharp contrast 
with standard broad transitions, the frequency Doppler 
broadening kvz remains larger than F. As a direct con- 
sequence, the optical depth b depends on the exact longi- 
tudinal velocity distribution f{vz) (see Eqs. ^ and (|3])) 
which are not necessarily gaussian [l^ . Since we have 
not access to f{vz) but only to Vz or T, one has enough 
control to assert the 6 <C 1 limit, thus the occurrence of 
the self-gravity regime, but we can perform only qualita- 
tive tests of our theory described above. 

At t = —50ms, the MOT cooling laser beams are 
turned off, leaving the trapped atomic cloud in an out- 
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FIG. 1: Upper part : Typical temporal evolutions of Lz the 
rms longitudinal size of the atomic cloud for three different 
fD beam intensities, 7 = (blue circles), I = 0.0157s (green 
circles) and 7 = 0.17s (red circles). The laser detuning is 
5 = — 5r for all curves. Lower part : The center of mass 
position of the atomic gas without the ID lasers (7 = 0). The 
y axis origin is arbitrary. 



of-equilibrium macroscopic state. Without the ID lasers, 
we observed a weakly damped oscillation of the breath- 
ing mode and of the center of mass position (blue circles 
in Fig. [T]). One notes that damping is caused by anhar- 
monicity of the dipole trap and not by thermalization of 
the gas which is negligible on the experimental timescale. 
In presence of the ID laser beams, overdamped or un- 
derdampcd oscillations of the cloud are observed. Both 
regimes will be now analyzed in detail. 

Let us first consider the stationary state in the over- 
damped situation (red circles in Fig. [1]). After the tran- 
sient phase {t < 30 ms), the rms longitudinal size of 
the atomic gas reaches a plateau at a minimal value of 
~ 120 fim with T ~ 2 fiK. The slow increase of the 
cloud's size after the plateau {t > 150 ms) goes with an 
increase of the temperature up to 4 /xK at the end of the 
time sequence. The origins of the long time scale evolu- 
tion are not clearly identified, but it is most likely due 
to coupling of the longitudinal axis with the uncooled 
transverse dimensions because of imperfect alignment of 
the ID laser beams with the longitudinal axis of the trap 
and nonlinearities of the trapping forces. At the plateau 
where temperature is around 2 /xK, we stress that, ac- 
cording to Eq. ([8]), the non- interacting gas is expected 
to have a rms longitudinal size of = L„i ~ 370 fim. 
Hence, here we clearly observe a compression of the gas, 
by a factor of three, due to the attractive interaction 
induced by the absorption of the ID laser beams. More- 
over, the estimated optical depth is 5 ~ 0.4. We then 
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FIG. 2: Density linear distribution in self-gravitating limit. 
The black dots are the experimental data. The profiles were 
symmetrized to improve the signal to noise ratio. The blue 
dotted and red dashed curves are fits of the experimental data 
with respectively, a gaussian and a sech^ profile. The fits are 
performed taking the distribution width and height as the two 
free fitting parameters. 



almost fulfill the two previously mentioned conditions 
— 6^1 and ^ Lni — for being in the ID self- 
gravitating regime. We note that the corresponding ex- 
perimental linear density distribution (black dots in Fig. 
[2]) has longer tails than the standard Gaussian profile 
(blue doted curve in Fig. [5]). Moreover, we find an ex- 
cellent agreement with the expected sech^ profile (red 
dashed curve in Fig. [2]). We have also checked that 
in absence of the ID laser beams, the experimental lin- 
ear density distribution corresponds to the expected pro- 
file of a non interacting gas in our dipole trap having a 
Zfi = 1.2(1) mm Raylcigh length. 

In the self-gravitating regime a 1/7V dependency of 
Lz is expected at fixed temperature (see Eqs. ^ and 
dSl)). Fig. [3] shows that the cloud's size Lz is in agree- 
ment with this prediction for two temperature ranges: 
1.5(2) /iK (blue circle) and 2.1(2) /iK (red star). Linear 
fits correspond to the blue dashed line for 1.5(2) /iK and 
the red dashed line for 2.1(2) /iK. The fitting expression 
is iV^ = 0.1 /L^ + 0,2, where ai and 02 are free parame- 
ters. The presence of the holding trap is revealed when 
N goes to zero by the finite value of L^- The fit gives 
Lz = \J~o.\la2 — 400 fj.m ioi N = 0, slightly larger but 
still in reasonable agreement with the expected value of 
Lni at these temperatures. The 1/A^ dependency of Lz 
in the self-gravitational regime is a clear signature of the 
long range nature of the interaction. 

Let us now consider the evolution of the trapped cold 
cloud in the underdamped situation (as an example see 
green circles in Fig. [T]). Without the ID lasers, the ratio 
of the eigenfrequencies of the breathing mode Wbr and 
the center of mass uJz is found to be close to two, as 
expected for a non-interacting gas in a harmonic trap. 
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FIG. 3: Dependency of the longitudinal size of the cloud with 
the number of atoms for 5 — 5.7(5)r and I — 0.31s- The 
blue circle (red star) data points correspond to temperature 
1.5(2) /iK (2.1(2) ^lK). The optical depth is in the range of 
0.6-0.2 according to atoms number variations. The blue and 
the red dashed lines are linear fits. 
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small optical depth, we check, using the Fokker-Planck 
equation with the ansatz (O, that Eq. (fTO|) with a = 
still provide a reasonable approximation for the breath- 
ing frequency. Fig. |4] summarizes the comparisons be- 
tween the measured ratio (wbr/<^z)^ and the predictions 
deduced from the relation pil)l . p is computed from the 
experimental data in the stationary state. We expect 
a = 0, however to judge the nature of the long-range at- 
tractive interaction, three plots respectively for a = 0, 1 
and 2 are shown. If the a = 2 case can be excluded, 
the experimental uncertainties does not allow to clearly 
discriminate between a ~ and a = I. 

To conclude, we give some experimental evidences of 
an ID gravitational-like interaction in a Strontium cold 
gas, induced by quasi resonant contra-propagating laser 
beams. First, we show that in the self-gravitating limit, 
the density distribution follows the expected sech^ pro- 
file. Moreover, the scaling of the cloud size with the 
number of atoms follows the predicted 1/N law. Finally, 
the modification of breathing frequency of the cloud, due 
to the long range interaction, is correctly described by 
a self-gravitating model. This suggests interesting con- 
sequences: by contrast with the ID case, a 2D self- 
gravitating fluid undergoes a collapse at low enough tem- 
perature, or strong enough interaction. Hence, it is con- 
ceivable that an experiment similar to the one presented 
in this paper, in a pancake geometry, would show such a 
collapse. 
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FIG. 4: Comparison for a = 0, 1 and 2 of the experimental 
ratio (cjbr/t*Jz)^ and the predictions deduced from the relation 
(|10|l . The values of p are measured on the experiment. 



As an example the blue curve, shown in Fig. (TJ gives 
Wbr/'^2 = 1.9(1). If now the attractive long-range interac- 
tion is turned on, some modifications of Wbr are expected 
whereas uj^ should remain unchanged. More precisely, 
assuming a power law two-body interaction force in the 
gas l/r", it is possible to find a simple relation for Wbr 
in the weak damping limit [l^ : 



cjbr=W2((3-a)(p-l)+4) = 



(10) 



This formula relates Wbr to a, the index of the interaction 
potential, and p = L'^jLl. Eq. (fTO)) is derived in [13] as- 
suming a velocity independent interaction term, which 
would be obtained by linearizing the radiation pressure 
force H]) in velocity. Since 0.3 < kvz/\S\ < 1.2 for our 
experimental data points shown in Fig. 21 the force can- 
not reasonably be linearized. However in the limit of 
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